What is claimed is: 



1 . A method for inde}|cing feature vector data space 
comprising the step of: 

(a) adaptively approximating feature vectors on the basis of 
statistical distribution of feature vector data in the feature vector data 
space. 

2. The method of /claim 1 , wherein the step (a) further 
comprises the steps of: 

(a-1) measuring the statistical distribution of the feature vector 
data in the feature vector data space; 

(a-2) estimating marginal distribution of the feature vector data 
using the statistical distribution; 

(a-3) dividing the estimated marginal distribution into a plurality 
of grids in which a probabi ity of disposing the feature vector data in 
each grid is uniform; and 

(a-4) indexing the feature vector data space using the divided 

grids. 

3. The method of claim 2, further comprising prior to step (a- 
4), the step of updating the grids on the basis of a previous probability 
distribution function and an updated probability distribution function, 
when new data is received. 
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4. The method of 
comprises indexing using vector 



claim 2. wherein step (a-4) further 
approximation (VA) files. 



5. The method of claim 2, wherein a number of the plurality 
of grids is determined by a number of bits assigned to the dimension. 

6. The method of claim 2, wherein step (a-2) further 
comprises the steps of: 

(a-2-1) defining a probability distribution function using a 
weighted sum of a predetermilied distribution function; and 

(a-2-2) obtaining an estimated probability distribution function by 
estimating predetermined parameters using the probability distribution 
function defined in the step (a-2-1). 



7. 



The method of 



claim 6. wherein step (a-2-2) further 
comprises obtaining the estimated probability distribution function by 
estimating the predetermined parameters using all N predetermined 
data in each estimation, wherein N is a positive integer, on the basis of 
an expectation-maximization algorithm using the probability distribution 

(a.2-1). 



function defined in the step 



8. 



The method o 



f claim 6. wherein the predetermined 



distribution function is the Gaussian function 
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9. The method of claim 6, wherein the probability distribution 
function of step (a-2-1) is a one-dimensional signal, p(x), wherein 

N 

p(x) = ^p(x\j)PU) , and wherein p(x\j) is defined as 



p(x 



2^ 



5 wherein coefficient P(j) is a mixing/ parameter that satisfies the 
following criterion 0<P(j)<l and ^P(j) = l. 

f 

10. The method of claim 6, wherein the estimated probability 
distribution function of step {a-Jr2) is obtained by finding <I)y,y=1,...,A/f, 



which maximizes <D(<D,,...,<Dj^)/= flP(v[/]|(a),,...,0^)) , where 
parameters y[f\, /=1,....A/, Is a given data set. 



1 1 . The method of claim 10. wherein the estimated 

I 

parameters of step (a-2-2) are updated according to the following 



equations 



i;^o|v[p'v[/] 



a'" = 



Epo1v[/])'(v[/]-//;.)^ 

= ^ h . and 
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representing a number of iterations. 



1 2. The method of claim 1 1 
parameter set of step (a-2-2) using N 



{pijY /^j .(o-^y}, and the updated parameter set for new data v[N+1], 



wherein t is a positive integer 



wherein the estimated 
data v[l] is given as 



coming in, is calculating using the fol 



//f^=//; + ^;^^(v[A^+,i]-//;). 



owing equations: 



-(POjv[7V + l])-P(7)^),and 



PO|v[Ar + l])' ^ 



1 3. The method of claim 



1 1 , wherein the step (a-2-2) further 



comprises: 

measuring a change of a probability distribution function which is 



defined as 



dx 



for each dimension, wherein a 
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previous probability distribution function is Poid(x), and an updated 



probability distribution function is Pne^,(x) \ and 



updating an approximation for the dimension if p is larger than a 
predetermined threshold value. 

14. The method of claim 2, wherein step (a-3) further 
comprises dividing a probabilitj distribution function into the plurality of 
grids to make areas covered by each grid equal, wherein the plurality 
of grids have boundary points/defined by c[l\, /= 0....,2*^, where b is a 
number of bits allocated and wherein the boundary points satisfy a 

criterion, pix)dx = — ifrl^^^ p(x)dx, and wherein the estimated 



probability distribution function is p{x) . 
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